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We present results of the electronic band structure, Fermi surface and electron transport properties 
calculations in the orthorhombic n- and p-type SnSe, applying the Korringa-Kohn-Rostoker method 
and the Boltzmann transport approach. The analysis accounted for the temperature effect on 
crystallographic parameters in Pnma structure as well as the phase transition to CmCm structure 
at T c ~ 807 K. Remarkable modifications of the conduction and valence bands were notified upon 
varying crystallographic parameters within the structure before T c , while the phase transition mostly 
leads to the jump in the band gap value. The diagonal components of the kinetic parameter 
tensors (velocity, effective mass) and resulting transport quantity tensors (electrical conductivity a, 
thermopower S and power factor PF) were computed in a wide range of temperature (15 — 900 K), 
and hole (p— type) and electron (n— type) concentrations (10 17 — 10 21 cm -3 ). SnSe is shown to have 
a strong anisotropy of the electron transport properties for both types of charge conductivity, as 
expected for the layered structure, with the generally heavier p-type effective masses, comparing 
to n-type ones. Interestingly, p-type SnSe has strongly non-parabolic dispersion relations, with 
the ’pudding-mold’-like shape of the highest valence band. The analysis of a, S and PF tensors 
indicates, that the inter-layer electron transport is beneficial for thermoelectric performance in 71 - 
type SnSe, while this direction is blocked in p-type SnSe, where in-plane transport is preferred. Our 
results predict, that the 77-type SnSe is potentially even better thermoelectric material than the p- 
type one. Theoretical results are compared with the single crystal p-SnSe measurements, and a good 
agreement is found below 600 K. The discrepancy between the computational and experimental data, 
appearing at higher temperatures, can be explained assuming an increase of the hole concentration 
vs. T, which is correlated with the experimental Hall data. 


I. INTRODUCTION 

Thermoelectric (TE) conversion in crystalline solids 
constantly attracts the interest of scientists, not only 
due to the increasing performance of energy harvesting 
systems, but also due to the fact that the conventional 
and ’well-known’ thermoelectric materials, such as PbX 
(.X =S, Se, Te)PEl or Bi 2 Te 3 (e.gP), still surprise by 
their novel and remarkable physical behaviors. It is com¬ 
monly accepted, that TE figure of merit zT = ctS 2 /kT 
well captures basic transport properties of the material 
at given temperature T, i.e. electrical conductivity (cr), 
Seebeck coefficient (S) and thermal conductivity (k), and 
is conveniently expressed in dimensionless units. As all 
these transport quantities apparently depend on tem¬ 
perature and carrier concentration (n or p), the maxi¬ 
mum of zT is expected when properly correlating both 
intrinsic electron transport properties of the system (e.g. 
the band gap magnitude) with temperature range, and 
doping level with the hole or electron concentration, to 
achieve p— or 77 —type materials, respectively. But, get¬ 
ting a better insight into atomic-level connections among 
crystal stability, electronic and lattice properties of TE 
systems generally allow for a more convinced interpreta¬ 
tion of the complex transport phenomena. 

Recent experimental work 1 reported zT ^ 2.6 at 
T ~ 920 K, along one of the axis in the single crystal 
p— type SnSe, which classified this well-known semicon¬ 
ductor as a very promising TE material. However, fur¬ 


ther result^ 8 1 revealed much lower zT for poly crystalline 
samples, showing that SnSe is a rather complex system. 

SnSe is a layered material, with an orthorhombic sym¬ 
metry of the unit cell, and its crystal structure strongly 
depends on temperature. First of all, SnSe exhibits a 
crystallographic phase transition at T c = 807 Kp Below 
T c , its unit cell can be described in the Pnma simple 
orthorhombic space group No. 62, with 4 chemical for¬ 
mulas (8 atoms) in the unit cell (see, Fig. [I]). Crystal 
structure parameters of this ’low temperature’ (i.e. be¬ 
low T c ) phase depend on temperature as well 6 . Since the 
temperature dependence of the transport properties (e.g. 
thermopower) is discussed in our work, those changes in 
crystal structure were also taken into account. Accord¬ 
ingly, three cases are considered: low-temperature (LT) 




FIG. 1: (Color online) Crystal structure of SnSe (Se - 
large, green balls, Sn- small, magenta). Left: LT phase 
(Pnma at 295 K); center: MT ’’phase” (Pnma at 790 
K); right: HT phase (Pbmm above 807 K). 
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phase, described by the experimental crystallographic 
data taken at T expt = 295 K, and used in our calcula¬ 
tions to represent SnSe in the temperature range T ca i c = 
10 — 550 K; mid-temperature (MT) phase (T expt = 790 K 
and T ca i c = 550 — 807 K); and high-temperature (HT) 
phase (T expt = 829 K and T ca i c > 807 K), as displayed 
in the Table [TJ Above T c , the high-temperature phase 
has an orthorhombic c-base centered CmCm space group 
(No. 63). The transition from the simple to the centered 
unit cell reduces the number of atoms in the primitive 
cell, from 8 to 4, changing the shape and size of the Bril- 
louin zone (BZ). Thus, it becomes impossible to directly 
compare the electronic dispersion relations between the 
LT/MT and HT phases, if using the HT centered unit 
cell. Therefore, to allow for a better understanding, and 
easier analysis of the role of the phase transition in the 
evolution of the electronic structure and transport prop¬ 
erties of SnSe, we transformed the HT SnSe centered unit 
cell into the equivalent Pbmm simple orthorhombic one 
(No. 51). Pbmm unit cell has the same number of atoms 
and similar alignment of the unit cell axes, as the struc¬ 
ture before the transition (i.e. a axis is the longest, b 
is the shortest, and c is in between, see, Fig. [I]). With 
this transformation, the Brillouin zone and location of 
the high symmetry points remain similar in both phases 
(small changes are only in BZ dimensions, due to lattice 
parameters’ variation). In our computations, the actual 
values of lattice parameters for HT phase were taken from 
the high temperature (T = 829 K) neutron measure¬ 
ments 6 , and are shown in the Table. [I] It is worth noting 
that the crystal structure of SnSe evolves smoothly with 
temperature (even while crossing T c ) and the unit cell di¬ 
mensions, as well as interatomic distances, change rather 
continuously with temperature (see, Table [I] and Ref. [6]) . 
The crystal structure parameter, which exhibits the most 
rapid change at T c , is the £ parameter positioning Sn and 


TABLE I: Crystallographic datsP of the SnSe 
compound, measured at different temperatures T expt , 
used in our calculations for selected temperature ranges 

Tcaic • 


Phase, 

T range 

Space 

group 

Lattice 
const. (A) 

Atomic positions 

Sn Se 

Low-T. (LT) 
Texpt = 295 K 
Tcaic < 550 K 

Pnma 

No. 62 

a = 11.501 

b — 4.153 

c = 4.445 

x: 0.6186 0.3551 

y: 0.25 0.25 

z: 0.1032 0.4818 

Mid-T. (MT) 
Te X pt = 790 K 
T ca i c =550-807 K 

Pnma 

No. 62 

a = 11.620 

b = 4.282 

c = 4.334 

x: 0.6230 0.3553 

y: 0.25 0.25 

z: 0.0584 0.4878 

High-T. (HT) 
Texpt = 829 K 
Tcaic > 807 K 

Pbmm a 

No. 51 a 

a = 11.62 

b = 4.282 

c = 4.293 

x: 0.6248 0.3558 

y: 0.25 0.25 

z: 0.00 0.50 


Se atoms. 

Recently, first principles calculations of transport 
properties of p-type SnSe were reported in Ref. [9'. In our 
work we give more extended analysis of the electronic 
structure, effective masses and transport properties of 
both, n- and p-type, SnSe. The constant relaxation time 
approximation is used, and phonon drag effects are not 
discussed. 

The paper is organized as follows. In Sec. [IT] theoret¬ 


ical and computational details are presented. Sec. |III A 
describes in details the electronic band structure and the 
Fermi surface of the valence and conduction states, to¬ 
gether with the effectiv e mass and the transport func¬ 
tion analysis. Sec III B| discusses anisotropic (single crys¬ 
tal) and isotropic (polycrystalline-like) transport prop¬ 
erties in function of temperature and carrier concentra¬ 
tion, and comparison with experimental data is done. 
Sec. [TV] contains summary and conclusions of the work. 
In Appendix [A] formulas for averaging of the anisotropic 
transport coefficients, in order to compute their isotropic 
(polycristalline) analogs, are derived. Appendix [B] con¬ 
tains complementary data. 


II. THEORETICAL DETAILS 


A. Thermopower 


The Boltzmann transport theorjP^®^, successfully ap¬ 
plied to study transport properties of various crystalline 
materials 13 is used to calculate the energy-dependent 
electrical conductivity cr(<f) (so-called transport func¬ 
tion, TF): 


a(£) = e 2 y] f jLr ra (k)v ra (k) <g> v ra (k)<5(<? - <f„(k)). 

n J 

(i) 

Symbol 0 represents outer product (Kronecker product) 
of two vectors, v n (k) = V/c<f n (k) and r n (k) is a velocity 
and life time of electrons on band <^ n (k), respectively. 
The transport function tensor has to be reformulated 
into the form, that is more convenient for the numer¬ 
ical computation. With the use of the constant relax¬ 
ation time approximation (r n (k) = To), and after chang¬ 
ing 3D k-space integration to the 2D surface integration 
(<f n (k) —)► S n (£)) it takes the form: 


(T T {g) = T 0 e ~Y, J 

n S n (g) 


dS w{S n {g))®w{S n {g)) 

4tt 3 |v(S n W)| ' 1 j 


TF is directly related to the macroscopic transport 
coefficients, like thermopower, electronic part of ther¬ 
mal conductivity or electrical conductivity. Within this 
approach!^ the two basic transport tensors (electrical 
conductivity a e and thermopower S) can be expressed 
as 


e=^ { °\ S 


1 jfW 


“Transformed from CmCm (No. 61) 




eT ££(°) 


(3) 
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where 


= j dS 


(4) 


The value of the chemical potential fi c = p c (T, n^) de¬ 
pends on the temperature (T), carrier co ncent ration and 
the type of conductivity (n, p) (see, Sec. II C| ). 

In case of the anisotropic structure, like the orthorhom¬ 
bic SnSe, transport function tensor (and thermopower as 
well) has three linearly independent elements, and in this 
work, diagonal elements (i.e. S xx , S yy , S zz ) are calcu¬ 
lated. For the polycrystalline sample without texture, 
those three elements have to be averaged in a proper way, 
to obtain the macroscopically isotropic thermopower and 
power factor (PF). The relevant formulas are given here 
and derived in the Appendix [Aj 


^avg 


^xx&xx T 


Syy&yy T S ZZ (T Z 


+ &yy + &z 


( 5 ) 


PFavg — ^{S X x (T xx + SyyVyy T S zz < J zz)- (6) 

Since the power factor, in the constant relaxation time 
approximation, depends linearly on the relaxation time 
r (not calculated here), consequently PF divided by r is 
presented. 


B. DOS and effective mass 


Using Eq. [9] and [10| one may discuss the direction de¬ 
pendence of the effective mass. In case of an orthorhom¬ 
bic structure with orthogonal axes, [Mij] tensor is diag¬ 
onal, and components m* x , m* y , ra* 2 are computed. 

The isotropic band effective mass is determined by the 
geometrical mean 


m *so = ^m* xx m* yy m* zz . (11) 

Both ways of effective mass calculations (DOS effective 
mass from Eq. [8] and band effective mass from Eq. © 
give the same results only for parabolic bands. Differ¬ 
ence of those two results gives opportunity to estimate 
the importance of non-parabolicity of the electronic band 
structure. 

The aforementioned effective masses correspond to the 
T — OK temperature. The temperature effects can 
be taken into account via the Fermi-Dirac distribution 
function in a similar way, as in the thermopower cal¬ 
culation in Eq. [3j The actual number of ’active’ (con¬ 
ducting) electrons, in temperature T, can be represented 

as n act i ve = f dS > g(S > ) and the effective mass of 

active electrons, as a function of temperature, can be 
determined from 

d£m*(£)g(£) (--g 


m*(T,n d ) = 


/ 



In a similar way to the transport function, the density 
of states (DOS) is calculated: 

/ \^ f dS 1 / l-r\ 

5( } = ^ J 4^iv fc «g> n (k)r (} 

S n (<?) 

With this definition g(S) has the units of eV -1 m -3 , i.e. 
it includes the volume of the unit cell. 

The DOS function is closely connected to the DOS 
effective mass, m*. Here, energy dependent m* is calcu¬ 
lated using the formula 1 ^ 

TODOs(^) = m e m^ os (S’) = h 2 \J 7T 4 (jr(^’)g , (^’). (8) 

Alternatively, effective mass can be computed by in¬ 
tegrating the effective mass tensor, over the isoenergetic 
surfaces: 


= m e m*j(S > ) = J [Mij]dS / J dS, (9) 
sw(<?) ' £»(<?) 

where the effective mass tensor is defined as 10 


[Mij]- 1 


1 d 2 S 
h 2 dkidkj 


( 10 ) 


where rid is the ca rrier concentration, at which m* is cal¬ 
culated (see, Sec. IIC). Note, that this analysis requires 
an assumption that effective mass is well defined, which 
for strongly non-parabolic bands may not be valid. 


C. Doping and chemical potential 


SnSe is an intrinsic semiconductor, where the Fermi en¬ 
ergy lies inside the gap. As already mentioned, transport 
properties of the intrinsically (defect-doped) p-type single 
crystal, polycrystal 7 , and p-type Ag-doped polycrystaP 
samples were reported. To simulate the behavior of the 
system after doping, we use the rigid band modePP In 
this approach, additional number rid of electrons or holes 
is the control parameter, added to mimic n— type (posi¬ 
tive rid) or p—type (negative rid) behavior, and chemical 
potential g c = p c (T, rid) needed in Eq. © is calculated 
using the formula 


n + n d = 


[ d£g {£)— 

J 1 


exp 


( ^-/x c (T,n d ) \ 

V k B T J 


(13) 


n is the total number of valence electrons in the system 
(which is 10 per f.u. in case of SnSe) and the integral is 
taken from the bottom of the valence bands. 
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D. Band structure computational details 

Electronic band structure calculations were performed 
using the full potential Korringa-Kohn-Rostoker 16 21 
(KKR) method, within the scalar relativistic ap¬ 
proach 22 23 . Local density approximation (LDA) 
parametrization of Perdew and Wang 27 was employed. 
The self-consistent cycle was repeated until the difference 
between the input and output potentials was less than 
1 mRy in any point inside the unit cell. Isoenergetic sur¬ 
faces Sni^) were obtained with the use of marching cube 
algorithm 24 on a mesh consisting of 80 x 80 x 80 voxels. To 
improve visualization, vertex normaP' and vertex color 
techniques were also used. All energy dependent func¬ 
tions were calculated with a resolution of 2.5 meV. The 
transport function was additionally interpolated between 
energy mesh points, using spline functions, which allows 
to obtain converged results for temperatures above 10 K, 
in the concentration range of 10 17 — 10 21 cm -3 . 

III. RESULTS 
A. Electronic Structure 

Electronic dispersion relations are shown in Fig. [2} 
for the LT (Fig. [2p MT (Fig. |Jd) and HT (Fig. |2p) 
phased of SnSe. As expected, local density approxima- 
tion (LDA), used in this work, underestimated the band 
gaps. Since the gap value is the important parameter in 
transport properties calculations, especially in elevated 
temperatures, the computed band gaps were expanded 
to mimic the experimental ones (see Tab. [TT| ) . For the LT 
phase, the band gap value was set to E® xp = 0.86 e\®. 
In MT and HT phases, where experimental data are not 
available, the calculated gap values were extrapolated 
proportionally as in the LT phase, i.e. we assume that 
the LDA underestimation is proportional to the real gap 
value. It is worth noting, that the calculated values of 
the energy gaps for the LT and MT phases are almost the 
same, whereas the HT phase has much smaller gap (see, 
Tab. 0 - To verify whether the reduction of the band gap 
is related to the change in atomic positions or the unit 
cell parameters, additional calculations were performed 
for the MT phase, using (i) the MT atomic parameters 
and HT lattice constants, and (ii) MT lattice constants 
and HT atomic parameters. The resulting LDA band 
gaps were: 0.465 eV and 0.355 eV, respectively, thus we 
found, that the gap value is controlled mainly by the Sn 
and Se atomic positions. Recent calculations 9 ^, using the 
GW method, reported E g = 0.829 eV for the LT phase 
and 0.46 eV for HT, thus reduction of E g after the phase 
transition was also found. 

As we explained in the Introduction, due to the trans¬ 
formation of the HT crystal structure to the simple or¬ 
thorhombic Pbmm space group, Brillouin zones in all 
(LT, MT and HT) phases in our calculations have the 
same shape and high symmetry points (BZ is shown 


schematically in Fig. [2]). The reciprocal space x, y , 2 axes 
are parallel to the a, 5, c real space directions. Let us 
start with some general comments of the band structure, 
before going into the detailed analysis of the effective 
masses and the transport properties. In the LT phase, 
the highest valence band (VB) appears in the T — Z di¬ 
rection, parallel to the z axis. It is worth noting, that 
in the LT phase, this highest VB has a ’pudding mold’- 
like shape, reported to be beneficial for thermoelectric 
performance in other systems 28 ! For the higher temper¬ 
ature structures (MT and HT), valence bands are much 
different, which points out that electronic band structure 
is sensitive to temperature, as simulated by the unit cell 
changes. The ’pudding mold’-like VB changes the shape 
to become more parabolic, and the second VB at T-Y 
direction, aligns to band in T-Z direction. The lowest 
conduction band (CB), seen in the T — Y direction in 
all phases, is parabollic-like. Above the band gap we see, 
that the whole band structure in the T — Z direction con¬ 
siderably moves down, as going from LT to HT, which is 
correlated with the shortening of the corresponding unit 
cell c axis. In the LT structure, the in-plane unit cell pa¬ 
rameters, b = 4.153 A and c = 4.445 A are significantly 
different, which in the reciprocal space is reflected in a 
different alignment of the T — Y and T — Z bands. After 
the phase transition, as well as in the MT case, b ~ c, 
and the energy locations of CB’s minima and VB’s max¬ 
ima become similar between those two directions. What 
seems quite surprising, the differences in the band struc¬ 
tures between the LT and MT (both before the phase 
transition) are much larger than between the MT and 
HT structures (i.e. induced by the phase transition). 
This shows that the electronic structure evolves contin¬ 
uously with temperature, however the band gap changes 
abruptly at the phase transition (see, Tab. [TT|) , which can 
result in rapid changes in the transport properties. 

It is also interesting to analyze the bands in. Z — U 
and T — X directions, which represent the real-space a 
direction, i.e. contributing the charge transport between 
the SnSe layers. In all the cases, for the VB part, the re¬ 
sulting bands are very flat, which is quite intuitive, since 
we expect hampered charge propagation in this direc¬ 
tion. On the other hand, conduction band in T — X is 
very steep and almost linear. All these band structure 
features are reflected in the transport properties of SnSe, 
as discussed below. 


TABLE II: Calculated (LDA) and experimental 1 values 
of the energy band gap (in eV) in SnSe for the three 
structures, LT, MT and HT. 


Structure 

LDA 

Experimental 

Used 

LT 

0.474 

0.86 a 

0.86 

MT 

0.487 

no data 

0.87 

HT 

0.350 

no data 

0.64 
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FIG. 2: (Color online) Electronic band structure of the LT (a), MT (b), and HT (c) phases of SnSe, computed along 
the high symmetry directions, as shown in the orthorhombic Brillouin zone. Horizontal lines mark the Fermi energy 

positions (/i c at T = 0 K) 

for the electron/hole concentration (in cm -3 ) of n — 10 18 (red), n — 10 19 (blue), n — 10 2 ° (green) and n — 10 21 (orange), 

respectively. 








FIG. 3: (Color online) Fermi surfaces of the LT (a,d), MT (b,e) and HT (c,f) phases of SnSe, for the n —type (a-c) 
and p —type (d-f) doping, for the carrier concentration (in cm -3 ) 10 19 (left sub-panel) and 10 2 ° (right sub-panel). 

Electron velocities (in m/s) are represented by a color scale. 


1. n-type doping 

For the n-type doping (whatever the LT, MT and 
HT structures) at electron concentrations < 10 19 cm -3 , 
Fermi level reaches one ellipsoidal electron pocket (with 
two-fold degeneracy) in Y-Y direction (see Fig. [2] and 
Fig. [3] a,b,c). The band, forming this pocket, has a reg¬ 
ular, parabolic-like shape, and isotropic effective mass 
equals m* so ~ 0.3 — 0.4 in LT, and m* so ~ 0.2 — 0.3 in 
MT and HT case, respectively. The highest values of m* 
are found along the x-direction (see Fig. ,d,f), where 
they quickly exceed m* x ~ 0.5. Therefore, velocity inte¬ 
grated over the Fermi surface (i.e. transport function), 


is 3 times lower in the x-direction than in the y- and 
^-directions (see. Fig. j5j),d,f), and low x-direction elec¬ 
trical conductivity is expected (since a e ~ cr (&?)). This 
behavior corresponds to our expectations of the lower 
electrical conductivity perpendicular to SnSe layers (in¬ 
direction is along a-axis). Electron transport properties 
change dramatically when the Fermi level reaches fur¬ 
ther electron pockets; one in HT and MT at the Y point, 
at ^ 10 19 cm -3 , and five in the LT structure (one cen¬ 
tered at the T point and four in the T-T direction) at 
~ 10 2 ° cm -3 . At these concentrations, the aforemen¬ 
tioned steep linear band in Y — X is activated. In con¬ 
trast to the previous case, these new pockets have the 
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FIG. 4: (Color online) Effective mass of the LT (a,b), 
MT (c,d) and HT (e,f) phases of SnSe, calculated using 
Eq.|8j(m D os), Eq. [9] (m xx , m yy and m zz ) and Eq. 
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{miso)- In each case, zero energy is fixed to the middle 
of the energy gap. Black vertical lines show the Fermi 
level for various carrier concentrations (in cm -3 ) 

1 x 10 18 , 1 x 10 19 and 1 x 10 2 °. Left and right columns 
present results for the p-type and n-type doping, 
respectively. 


FIG. 5: (Color online) Transport Function (TF) of 
valence and conduction bands in LT, MT and HT 
phases of SnSe. Zero energy is fixed in the middle of 
energy gap, as in Fig. [2] Black vertical lines show Fermi 
level for various carrier concentration (in cm -3 ): 

1 x 10 18 , 1 x 10 19 and 1 x 10 2 °. Left and right columns 
present results for p-type and n-type doping, 
respectively. 


smallest m xx (see Fig. [4|, while the largest mass tensor 
component is now m zz . In consequence, the average ef¬ 
fective mass, represented by m* so or m ^ 0s , raises rapidly 
to 1.0 — 1.5 (for n ~ 10 2 ° cm -3 in LT ), and 0.6 — 0.9 
(for n ~ 10 19 cm -3 in HT and MT). Thus, the transport 
function tensor component a xx (Fig. changes from 

the smallest, in energy range corresponding to the elec¬ 
tron concentration below 10 19 cm -3 , to the largest, at 
energies corresponding to 10 2 ° cm -3 . The rapid increase 
of the g xx (connected to the alignment of the five pockets 
at the same energy) should favor the high thermopower 
S XXl since S ~ d(lna)/dS > . In case of the n-type doping, 
in the LT, MT and HT phases, after reaching certain 
concentration between 10 19 cm -3 (HT) - 10 2 ° cm -3 (LT, 
see, Fig. [5]) two types of electrons are involved in the elec¬ 
trical transport along the x-direction: first, those having 
low mass and high velocity, providing high conductivity, 
and second, heavy electrons that are needed to achieve 
high thermopower. In such a case large power factor can 
be expected. 


DOS effective mass, plotted as a function of the car¬ 
rier concentration and temperature (temperature effects 
according to Eq. 12), is shown in Fig. [6j For a sake of 
completeness, we present results for all the phases (LT, 
MT and HT) in the wide 15 — 900 K temperature range. 
Appearance of the additional electron pockets, that is 
manifested in the rapid rise of the effective mass, have 
sharply determined concentrations only at very low tem¬ 
peratures. For T = 300 K, due to temperature blurring 
of the Fermi-Dirac statistics, m* changes gradually with 
a pronounced bump above 10 2 ° cm -3 in the LT struc¬ 
ture. At temperature c.a. T = 800 K or higher, elec¬ 
trons, which actively conduct, presumably come from the 
deeper lying bands, that have strongly non-parabolic dis¬ 
persion relations. Thus, effective mass is not well defined, 
and the results must be treated with caution. Generally, 
in all the phases, DOS effective masses rise with the con¬ 
centration, from approximately 0.3 m e at n = 10 18 cm -3 
to 2 — 3 m e at n = 10 21 cm -3 . 
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2. p-type doping 

In the p-type LT-structure, up to 10 20 cm -3 , the Fermi 
level penetrates hole pockets located in Y -Z direction 
(Fig. §1). These pockets have non-ellipsoidal shape, 
which is caused by the ’pudding-mold’-like band, hav¬ 
ing two maxima (see Fig. [2^i). In this case, effective 
mass cannot be well defined (see Fig. [ 6 ^ 1 ) , since such a 
band exhibits both, electron-like, and hole-like proper¬ 
ties, due to the convex and concave curvatures. However, 
the Fermi surface shape indicates, that the x-direction ef¬ 
fective mass possesses the highest value (elongated shape 
in ^-direction), which, in real space, corresponds to the 
a-axis direction (perpendicular to the SnSe layers). This 
behavior obviously affects the carrier velocity, integrated 
over the Fermi surface (see, Fig. [5^). The component a xx , 
similar to n-type case, is again the lowest and the highest 
values of TF are detected in z- (close to VB edge) and 
p-directions (well below VB edge), both being parallel to 
the SnSe atomic layers. Interestingly, at higher temper¬ 
atures (MT and HT phases) the bands in Y-Z direction 
do not have inflection, and the effective mass is therefore 
well defined. At higher concentrations ( 10 20 — 10 21 cm -3 ) 
in LT , the Fermi level also reaches holes from the Y-Y 
direction (see Fig. [2|i and Fig. |3^i) . In HT (and also MT), 
the axial anisotropy of the transport properties of SnSe 
(rn* and TF) is nicely visible: very large effective mass 



FIG. 6 : (Color online) DOS effective mass in function of 
the carrier concentration and temperature for the LT, 
MT and HT phases of SnSe. Left and right columns 
correspond to the p-type and n-type doping, 
respectively. 


and the smallest transport function is seen in ^-direction, 
compared to much more ’conducting’ y- and ^-directions 
(in-plane transport). The electrical conductivity cr e , xx is 
expected to be about three times lower (see Fig. [ 5 J 3 ) than 
in the other directions. 

Figs. [6] (left column) show p-type DOS effective mass 
in function of the carrier concentration and temperature. 
The LT structure, already at n = 10 18 cm -3 , has con¬ 
siderably large value of the ’bare’ mass (i.e. not affected 
by temperature blurring, T = 15 K curve), m* ~ 1.0, 
which is raising significantly above n = 2 x 10 19 cm -3 . 
At T = 300 K, similarly to the n-type doping, critical 
concentration, where additional bands start to have in¬ 
fluence on rn*, is blurred, and rn* is almost constant 
(rn* ~ 1.5) up to 10 2 ° cm -3 , then rising to 3 m e at 
n = 10 21 cm -3 . At higher temperatures, the deep and 
heavy valence bands contribute to the effective mass even 
at the lowest concentrations, however, we have to keep 
in mind, that integral in EqUU covers carriers being far 
from the gap, where bands are generally not parabolic, 
thus the characterization of bands in terms of the effec¬ 
tive mass may become inaccurate. 

Closing the effective mass discussion, for all the LT,MT 
and HT cases, n-type effective masses are smaller than 
P- type. 


B. Transport properties 

1. Thermopower 

Thermopower as a function of electrons/holes concen¬ 
tration for the three different temperatures are presented 
in Fig. [ 7 ] where data for the different temperatures were 



FIG. 7: (Color online) Thermopower, in function of the 
carrier concentration, for the p- (top panel) and n-type 
(bottom panel) SnSe, for T — 70 K, 300 K and 807 K. 
The electronic band structures of different phases were 
used to calculate thermopower: LT at 70 K and 300 K, 
MT and HT at 807 K. 


























































































































































obtained for the corresponding crystal structures. The 
complete results for all the phases at all the temperatures 
are shown in the Appendix [B] The most interesting fea¬ 
tures are seen for the n-type LT case, that exhibits strong 
anisotropy of the Seebeck coefficient, visible at T = 70 
and 300 K around n ~ 10 20 cm -3 . The n-type S xx tensor 
element has the largest absolute value; at low tempera¬ 
tures S xx starts to dominate at 1 x 10 2 ° cm -3 , where the 
Fermi energy reaches 5 electron pockets with high v x (see, 
red line in Fig^). At higher temperatures (i.e. for MT 
and HT), around n ~ 2 x 10 19 cm -3 , Ep reaches only one 
additional electron pocket (see, Fig. iu red line), and 
the advantage of reaching this pocket spreads over wider 
concentrations range (see, n-type S xx at 807 K in Fig. [7]). 
This is reflected as ^ 20% rise of S xx thermopower in the 
whole concentrations range. 

The anisotropy of p-type thermopower is less signifi¬ 
cant. Similar enhancement of the thermopower, as ob¬ 
served for the n-type, caused by the appearance of addi¬ 
tional bands near the Fermi level, is seen in p-type S yy 
component at p ~ 3 x 10 19 cm -3 . Two, non-parabolic 
pockets, have a high velocity of electrons in ^-direction, 
increasing rapidly at p ~ 3 x 10 19 cm -3 , which provides 
the increase of a yy derivative (see, green line in Fig. [5^). 
Significance of this effect is smaller, than in the n-type 
doping, and seen only in the LT phase, but it is sufficient 
to provide an increase of S yy over S xx and S zz (best seen 
at 300 K). 

As far as the bipolar effects are concerned, we can ob¬ 
serve, that p-type S xx element shows the strongest reduc¬ 
tion at the high temperature, which is caused by the large 
n-type S xx and a xx . The major detrimental effect at high 
temperatures is related to the reduction of the band gap 
during the phase transition (LT and MT vs HT, see Ta- 
ble[!!|. To visualize this influence of the gap reduction on 
the bipolar effect, we have plotted S at the same temper¬ 
ature of 807 K (phase transition temperature) computed 
for the MT and HT cases. Much smaller band gap in the 
HT phase leads to the much stronger bipolar effect, and 
the thermopower starts to decrease already for p < 10 19 
cm -3 and n < 2 x 10 18 cm -3 . Note also, that the bipolar 
effect critically depends on the value of the band gap, 
which here for the HT phase was adjusted in an approx¬ 
imate way, since we are not aware of the experimental 
value. 

Isotropic thermopower, calculated using Eq. [5| is 
shown as a color map in Fig. [8j and S(T) curves for se¬ 
lected carrier concentrations are collected in Appendix \B\ 
Fig. [l6j To allow for a comparison between the different 
phases, all three cases for complete temperature ranges 
are displayed. After averaging over three directions, in 
all the three phases, p-type thermopower is generally 
larger than n-type one. For n-type LT map, around 
n ~ 10 2 ° cm—3, and below 400 K, we observe the ab¬ 
normal bending of S', which increases with the carrier 
concentration. This is due to the rise in S xx element, 
discussed before. For the MT and HT phases, such effect 
is not observed, and the variation of S with the carrier 
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FIG. 8: (Color online) Isotropic thermopower for LT, 
MT, and HT phases of SnSe (see Tab. |T| for n- and 
p-type doping. 


concentration is monotonic (S decreases with nor p). At 
the highest temperatures and the lowest carrier concen¬ 
trations, we again observe a drop of S due to the bipolar 
effects, the strongest for the p-type HT phase. 


2. Power factor 

The power factor (PF= S 2 /cr), similarly to the elec¬ 
trical conductivity, cannot be directly calculated if the 
electronic relaxation time r is not known, however the 
discussion of PF/r can be still very useful in optimiz¬ 
ing the carrier concentration in the SnSe system. For 
completeness of the analysis, we show results for all the 
phases (LT, MT and HT) at T = 450 and 807 K, although 
one has to bear in mind that these phases describe SnSe 
in limited temperature ranges. 

Results of the calculations are shown in Fig. [9| The 
most striking observation is, that for the p-type SnSe, xx 
tensor element of the power factor (PF^) is the smallest 
one, in contrast to the n-type, where PF^ is the largest. 
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FIG. 9: (Color online) Power factor of the LT, MT and 
HT phases for p- (left column) and n-type (right 
column) SnSe in function of concentration for two 
temperatures (450 K and 807 K), calculated along three 
directions. 


The reason is that in the p-type, ^-direction inter-layer 
conducting channel is blocked due to the smallest TF 
(see, Fig. [5j red line), whereas in the n-type, it is acti¬ 
vated by the electron pockets with high v x . The smallest 
p-type PF xx , comparing to y- and ^-directions, remains 
in agreement with the experimental findings 1 1 

In the p-type LT phase, value of the thermopower is 
approximately isotropic (except for S xx at low concen¬ 
trations), thus overall anisotropy in the PF is due to the 
electrical conductivity: transport function element CFyy/r 
has the biggest value (in p > 5 x 10 19 cm -3 ) which makes 
PF the largest in the p-direction. For the high temper¬ 
ature structures, namely MT and HT, the power factor 
for y- and ^-directions becomes similar, and is the largest 
for the HT phase. In the n-type SnSe, PF/r curves look 
quite similar, while comparing the LT, MT and HT cases, 
with decreasing anisotropy (measured as FF xx /FF yy and 
FF XX /FF ZZ ) upon going from LT to HT. 

For isotropic (polycrystalline) material, the power fac¬ 
tor, calculated using Eq. [6j is mapped in Fig. [lOj The 
first observation is, that generally the n-type SnSe has 
larger values of PFs, than the p-type one, thus the heav¬ 
ily doped n-type SnSe can exhibit even better thermo¬ 
electric performance, and better polycrystalline zT over 
the broad temperature range. It seems that the exis¬ 
tence of the structural phase transition is beneficial for 
the power factor value only for p-type SnSe, where HT 
phase exhibits the largest values, but only for large car¬ 
rier concentrations, where the bipolar effects are not so 



FIG. 10: (Color online) Isotropic power factor for n- 
and p-type SnSe (LT, MT and HT phases) in function 
of the carrier concentration and temperature. All values 
in 10 12 WK- 2 m- x s- x . 


important. However, at higher concentrations, the zT 
value may not benefit from the phase transition, because 
of the increasing contribution of the electronic part of the 
thermal conductivity. 

The maps of PF also reveal that SnSe is best suited for 
TE application at temperatures T > 600 K and should 
be rather heavily doped, for both n- and p-type above 
10 2 ° cm- 3 . 

Note, that all the presented transport calculations re¬ 
sults and predictions are based on the constant relax¬ 
ation time approximation, and, in principle, electron- 
phonon or electron-impurity scattering may markedly af¬ 
fect them. Especially for the case, where both heavy, and 
light electrons are involved in electrical transport, scat¬ 
tering effects may be of increased importance, due to the 
different scattering rates for these two electronic bands. 
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C. Comparison with experiment 


temperature dependence to our extracted p(T) function. 


Experimental data were taken from Ref. |l] where the 
transport properties of the single crystal p-type SnSe 
were reported. Since our calculations were performed 
within the constant relaxation time approximation, only 
thermopower can be directly compared. Fig. [Tl] shows 
the measured Seebeck coefficient along a, 6, and c crys¬ 
tallographic directions 30 compared with two sets of cor¬ 
responding theoretical values, for hole concentrations 
p = 3 x 10 17 and 7 x 10 17 cm -3 . These carrier con¬ 
centrations were selected for comparison basing on the 
measured Hall coefficient of SnSe at T = 300 K, de¬ 
picted also in Fig. 12 r. Nevertheless, one has to bear 
in mind that due to the anisotropic crystal structure and 
Fermi surface, as well as the strongly non-parabolic band 
structure, Hall coefficient data do not reflect directly the 
carrier concentration. In Fig. 11 we see, that the exper¬ 
imental Seebeck coefficient well matches the theoretical 
values below 600 K, for p = 7 x 10 17 cm -3 , and then is 
underestimated in the mid-temperature range 700 K). 
The most intriguing experimental behavior is noticed for 
the thermopower, which remains constant with temper¬ 
ature above T c . For such high temperatures, and car¬ 
rier concentration around p ~ 10 17 cm -3 , bipolar effect 
should strongly decrease 5, even making S negative, as 
seen from our computations, shown in Fig.[7]and Fig. [16) 
Similar observations, with even faster decreasing S due 
to the smaller band gap, were reported from calculations 
in Ref. |9l 

One can attempt to explain such an uncommon situ¬ 
ation assuming the gradual increase of the carrier con¬ 
centration with temperature, since it could explain both, 
the decrease in S above 600 K, and saturation of S above 
800 K, with almost invisible bipolar effect. Experimental 
Hall measurements also show, that 1 /Rr rises almost 
two order of magnitudes (see, Fig. fl2k), between 600 
and 800 K, which supports the above-mentioned hypoth¬ 
esis. Moreover, for the high carrier concentration, the 
thermopower changes smoothly while crossing the phase 
transition temperature. Fig. [ 7 ] shows that only concen¬ 
trations larger than 1 x 10 19 cm -3 provide similar Seebeck 
coefficient at T = 807 K for MT and HT phases. Finally, 
measured electrical conductivity (see, Ref. |Tj) shows non¬ 
monotonic behavior with increasing temperature, i.e. it 
decreases from 300 K to 550 K, increases above 550 K 
and becomes constant above T c . Such a trend remains in 
line with the expectation of increasing carrier concentra¬ 
tion above 550 K, where the experimental and theoretical 
thermopower starts to deviate. 

To partly verify this hypothesis, we have plotted how 
carrier concentration should change with temperature, 
to reach the agreement between the experimental and 
calculated thermopowers! 31 As we can see in Fig. 12 3, 
the unusual, experimental behavior of thermopower, di¬ 
rectly results in the rising carrier concentration in SnSe, 
to more than 1 x 10 19 cm -3 around 800 K, and the mea¬ 
sured inverse Hall resistivity (Fig. Il2a), exhibits similar 



900 


FIG. 11: (Color online) Comparison of the calculated 
and experimental valued of the thermopower vs. 
temperature. Nominal hole concentrations for the 
theoretical curves are 3 x 10 17 (orange) and 
7 x 10 17 cm -3 (green). 



FIG. 12: (Color online) (a) Concentration derived from 
the measured Hall coefficient, according to 
tih = 1 /eRn)\ (b) nominal carrier concentration at 
which the calculated thermopower equals to the 
experimental value. All concentrations are given in 
cm -3 , experimental data are taken from Ref. |TJ Vertical 
line separates the temperature ranges, where theoretical 
thermopower of LT (below 550 K) and MT (above 
550 K) phases was used. 


IV. SUMMARY 

Electronic structure and thermoelectric properties of 
n and p-type SnSe were studied, using the KKR method 
and the Boltzmann transport approach, within the con¬ 
stant relaxation time approximation. We have shown, 
that the temperature evolution of the crystal structure 
within the Pnma phase (i.e. before reaching the struc¬ 
tural phase transition at 807 K), leads to the significant 
changes in the electronic band structure. On the other 
hand, the phase transition, occurring at 807 K, leads 
mainly to the abrupt change of the energy band gap 
value, whereas modifications of <f(k) curves are minor. 
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The effective masses, analyzed in function of the carrier 
concentration and temperature, indicate, that overall p- 
type masses are larger than n-type ones, and SnSe ex¬ 
hibits strong anisotropy of the electron transport prop¬ 
erties for both types of charge conductivity. 

P -type SnSe, computed for the room temperature crys¬ 
tal structure parameters, has strongly non-parabolic dis¬ 
persion relations, with the ’pudding-mold’-like shape of 
the highest valence band. Due to the flat bands, yield¬ 
ing large effective masses, the inter-layer electron trans¬ 
port seems to be blocked, resulting in small power factor 
and Seebeck coefficients, when comparing with the cor¬ 
responding values computed for the other directions (i.e. 
in-plane electron transport). 

The opposite situation was found in the n-type SnSe, 
where the inter-layer direction exhibits band struc¬ 
ture features, highly beneficial for thermoelectric perfor¬ 
mance. In this case, when critical concentration between 
10 19 cm -3 (HT) - 10 2 ° cm -3 (LT, see, Fig. [ 5 ]) is reached, 
or when high velocity electron band is thermally acti¬ 
vated at high temperatures, two types of carriers are 
present simultaneously in the system, i.e. more local¬ 
ized electrons, having large effective masses, and highly 
mobile electrons, possessing low effective masses and high 
velocity. This particular combination results in high ther¬ 
mopower and power factor in the inter-layer direction. In 
view of our results, the n-type SnSe may be even better 
thermoelectric material, than the p-type one. 

Our theoretical study confirmed the strong anisotropy 
of the electron transport in SnSe. This results in the 
much worse thermoelectric performance in the polycrys- 
talline material, in agreement with the experimental re- 
portd^. 

On the whole, the computational results, derived from 
the KKR method and Boltzmann equations, show quite 
good agreement with the measured data 1 below T ~ 
600 K. The experimentally observed changes in S', cr 
and Rh for T > 600 K suggest, that in the measured 
SnSe samples, the carrier concentration was gradually 
increasing, and the discrepancy between the theoretical 
and measured Seebeck coefficient, appearing at high tem¬ 
peratures (T > 600 K) can be corrected when the actual 
Hall concentration is taken into account in calculations. 
Nevertheless, the mechanism responsible for such a gen¬ 
eration of additional carriers, suggested in our discussion, 
remains unknown. 


Appendix A: Polycrystalline random grain model 


In a poly crystalline material, where grains have differ¬ 
ent crystalline orientation and are randomly distributed, 
the effective Seebeck coefficient, as well as the power fac¬ 
tor, are not simple averages of the single grain properties. 
To obtain the correct isotropic values, appropriate circuit 
model for thermoelectric generator is build (see Fig. 13). 

In the polycrystalline sample it is assumed, that the 
material consists of three types of small grains, with 
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FIG. 13: Circuit model of lumped (with R L ) 
thermoelectric material. 


thermopower Sk {k = 1,2,3) each, and electrical con¬ 
ductivities <Jk (see Fig. 14), distributed on a regular grid. 
This model can be simplified to the thermoelectric mod¬ 
ule made from three solid materials connected in parallel. 
On this basis, equivalent circuit model was made. 



FIG. 14: Model of the thermoelectric material, which is 
build up from three types of small, randomly 
distributed grains, with different transport properties. 


The thermoelectric module with one the type of mate¬ 
rial, having the thermopower S and electrical conductiv¬ 
ity <7, can be replaced with the electromotive force (emf) 
£ = SAT and r = crl/A, where l and A are length, and 
cross section of the module, respectively (see Fig. [T3] ) . 
Effective thermopower can be defined as voltage on the 
module, divided by the temperature difference AT, when 
Rl goes to infinity. 


~ _ U _ iRl _ £ Rl r l —^00 £ _ „ 

eff ~ AT ~ ~AT “ r + R l AT > AT ~ 

(Al) 

where i is the current in circuit. 

The power factor can be defined as a capability of en¬ 
ergy production of material with the cross section A, and 
length /, at the given temperature difference AT 


PFeff AAT 2 ’ 


(A2) 


where P is the power output of the source, when Rl goes 
to zero (short circuit power). 


P = £i = 


£ 2 

r + Rl 


Rl^O, £ 2 
r 


(A3) 


and therefore 


PF eff = 


PI 


AAT 2 



(A4) 
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In case of the equivalent polycrystalline circuit model, 
the effective thermopower is: 


_ U _ iRl 

eff ~ ~ 


(A5) 


From circuit theory 


^k 


Sk — iRl 

Tk 


and i = ik 


(A6) 


where k is 1, 2, 3. Now i can be found 


Zfc £k/ r k 
R L J2k l / r k + 1 


(A7) 


_ iRL _ Zlfc ^k/l’k _P R ]_, —^OO 

eff= AT = R l £, l/r k + 1 IRl A T > 

E/cV r fe AT E/c^/c 


(AS) 


Power of the source is now a sum of powers of all the 
sources 


P = £ k i 

k 


£k ~ iSkRL 

Tk 





(A9) 


and therefore 

(A10) 

where Ak = 1/3A (grains are randomly distributed). 


a) 





FIG. 15: (Color online) Thermopower LT (a), MT (b) 
and HT (c) phases at three different temperature and 
crystallographic direction. 


Appendix B: Extended data 


In this Appendix, additional figures are included, 
showing the complete data for the anisotropic (Fig. [I5| ) 
and isotropic (Fig. 16) thermopowers of the three consid¬ 
ered SnSe structures. 
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FIG. 16: (Color online) Isotropic thermopower in func¬ 
tion of temperature for LT (a), MT (b) and HT (c) for 
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concentrations in cm -3 . 
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